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ABSTRACT 

We confirm an eighth gravitational lens system in the CAmbridge Sloan Survey Of Wide ARcs 
in the skY (CASSOWARY) catalogue. Exploratory observations with the X-shooter spec- 
trograph on the Very Large Telescope (VLT) of the European Southern Observatory (ESO) 
show the system CS WA 5 to consist of at least three images of a blue star-forming galaxy 
at z = 1.0686, lensed by an apparent foreground group of red galaxies one of which is at 
z = 0.3877. The lensed galaxy exhibits a rich spectrum with broad interstellar absorption 
lines and a wealth of nebular emission lines. Preliminary analysis of these features shows the 
galaxy to be young, with an age of ^ 25-50 Myr With a star-formation rate of ^ 20 Mq yr~^, 
the galaxy has already assembled a stellar mass il/* ^ 3 x 10^ Mq and reached half-solar 
metallicity. Its blue spectral energy distribution and Balmer line ratios suggest negligible in- 
ternal dust extinction. A more in-depth analysis of the properties of this system is currently 
hampered by the lack of a viable lensing model. However, it is already clear that CSWA 5 
shares many of its physical characteristics with the general population of UV-selected galax- 
ies at redshifts z — 1-3, motivating further study of both the source and the foreground mass 
concentration responsible for the gravitational lensing. 

Key words: gravitational lensing: strong - galaxies: abundances - galaxies: evolution. 



1 INTRODUCTION 

Galaxies at redshifts z > 1 are generally too faint for their rest- 
frame ultraviolet (UV) spectra to be recorded with the high res- 
olution and signal-to-noise (S/N) ratios required for an in-depth 
study of their most important physical properties. Consequently, 
our knowledge of the high redshift galaxy population is still largely 
based on analyses of stacked spectra which average together obser- 
vations of many ga l axies in order to improve the S/N ratio (e.g. 
IShaplev et al] l2003l : IVanzella et al.l '2009 : We iner et"ZI 12009). A 
few exceptions are sources which, thanks to their fortuitous loca- 
tion behind foreground mass concentrations such as massive galax- 
ies, or groups and clusters of galaxies, are strongly lensed so that 
their fluxes are significantly boosted. Observations of what has be- 
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come the archetypal strongly lensed galaxy, MS 1512 -cB58 (cB58 
for short) at z = 2.7276 (P ettini et al. 2000; Teplitz et al.ll200(]| ; 
IPettini et al.ll2002l ; lBaker^alJi2004 ; iSiana et al. 2008), have pro- 
vided clear demonstrations of the power of this technique in giving 
insights into the interstellar medium and stellar populations of dis- 
tant galaxies with a level of detail which would otherwise be un- 
obtainable until the next generation of 30+ m optical-infrared tele- 
scopes. 

The la rge area of sky sa mpled by the Sloan Digital Sky Sur- 
vey (SDSS. lYorketal.ll2000h has allowed systematic searches fo r 
strongly lensed galaxies (e.g. lAllam et al.ll2007l ; iLin et al.ll2009h . 
with the dual interest of identifying and studying the most mas- 
sive galaxies in the lensing objects, and distant galaxies in the 
lensed sources. The Sloan Lens ACS Survey, or SLACSQ targets 
galaxies with emission lines at more than one redshift, or red- 
shifted nebular emission lines superimposed on the spectra of early 
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Figure 1. Left: 20 x 20 arcsec SDSS g, r, i colour composite image of the CSWA5 system. Right: r'-band image of CSWA5 recorded with a 20 s long 
exposure of the X-shooter acquisition and guiding camera, at position angle P.A. = 8 degrees. The scale of the right hand image is extended by a factor of two 
relative to the left one. In the right-hand figure, the lensed images of the source are labelled iml-im4, reflecting the sequence of light arrival times (on the 
assumption that im4 is one of the lensed images). Two foreground massive red galaxies are labelled LI and L2. Note that the blue image im4 is not clearly 
visible in the SDSS composite on the left. The X-shooter slit is overlaid on the right-hand image of the field; we have also indicated the apertures used in the 
extraction of the one-dimensional spectra of iml, im2, and im3, as described in Section|2] Coordinates of all six labelled objects are given in Table[T] 



Table 1. SDSS coordinates and magnitudes for the different components of the CSWA5 lens system. 



Object RA (J2000) 



Dec (J2000) 



LI 

L2 

iml 

im2 

im3 

im4 



12:44:51.35 
12:44:51.00 
12:44:51.22 
12:44:51.27 
12:44:51.23 
12:44:51.44 



01:06:43.0 
01:06:44.0 
01:06:40.1 
01:06:45.2 
01:06:44.2 
01:06:43.5 



22.34±0.50 
22.33±0.79 
22.07±0.20 
20.79±0.08 



20.61±0.04 
20.43±0.06 
21.74±0.06 
21.09±0.04 



19.12±0.02 
18.44±0.02 
22.48±0.17 
20.91±0.05 



18.44±0.02 
17.87±0.01 
23.85±0.67 
20.67±0.07 



18.02±0.04 
17.67±0.04 
21.75±0.45 
20.42±0.19 



Notes. The magnitudes quoted are SDSS model magnitudes and are considerably more uncertain than the formal errors quoted because 
of blending. No SDSS magnitudes are available for images im3 and im4. 

The coordinates of the objects were determined from the X-shooter acquisition image (see Figure[T), relative to the SDSS photometric 
source catalogue. Positional uncertainties are ~ 0.1 arcsec. 



type galaxies. These candidate gravitationally lensed systems are 
then followed up with high spatial resolution images from the Ad- 
van ced Camera for Surveys on the Hubb le Space Telescope (HST; 
e.g. lBolton et al.ll2006l : lwrilis et al.ll2005h . In a complementary ap- 
proach. The CAmbridge Sloan Survey Of Wide ARcs in the skY 
(CASSOWARY) targets multiple, blue companions around massive 
ellipticals in the SDSS photometric catalogue as likely candidates 
for w ide-separation gravitational lens systems (see lBelokurov et al] 
l2009l) . 

These and other search strategies have led to numer- 
ous recent discoveries of strongly lensed galaxies at z ^ 
1 which have b een followed up in detail ov e r a range of 
wavelengths (see Lemoine-Busserolle et al. 2003; Cabanac et al. 
2005; Belokurov et al. 2007; Coppin et al. 2007; Swinbank et al. 
2007; Cabanac etal. 20081 ; IStark et alj l2008i ; iFinkelstein et al. 
2009; Hainline et al. 200^; [Ouider et alj|2009l; [siana et al, .2009, ; 
Dessauges-Zavadsk v et al.l2010l ; Ouider et al. l2010l. and references 
therein). However, all of these studies have of necessity focused 
only on limited wavelength intervals, mostly in the rest-frame UV 
(redshifted to optical wavelengths) or the rest frame optical (red- 
shifted into the near-infrared). It is only with the recent advent of 
the X-shooter spectrograph on the Very Large Telescope (VLT) 
facility of the European Southern Observatory (ESO), that both 



wavelength ranges can now be recorded at once with good sensitiv- 
ity. The obvious advantages have been demonstrated in ou r recent 
study of CASSOWARY 20 (CSWA 20. IPettini et alj|2O10l) , a blue 
star-forming galaxy at Zcm ~ 1.433 lensed into a wide separation 
('^ 6 arcsec) Einstein Cross by a massive luminous red galaxy at 
Ziihs ~ 0.741 with velocity dispersion (Jions — 500 km s~^. 

The current CASSOWARY catalogu^fl includes seven spec- 
troscopically confirmed lensed galaxies at z > 1. In this paper we 
confirm an eighth case with X-shooter observations of CSWA 5. 
As can be seen from Figure [T] this multiple system consists of two 
foreground red galaxies one of which (L2) has an SDSS spectrum 
which shows it to be at Zabs ~ 0.3877, and four fainter sources, 
separated by ~ 1-5 arcsec. Three of these images, iml, im2 and 
im3, are shown by the observations reported here to be gravitation- 
ally lensed images of an emission line galaxy at Zcm ~ 1.0686, 
while the fourth image, not covered by our observations, remains 
to be confirmed spectroscopically. Coordinates and magnitudes of 
the different components of CSWA 5 are given in Table[T] However, 
we caution that the SDSS pipeline magnitudes, which are based on 
exponential profile fits, are likely to be significantly more uncertain 
than formal errors quoted, because of blending between the differ- 
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Figure 2. SDSS spectrum of galaxy L2 at z^bs = 0.3877 (shifted to the rest frame). The best fit stellar template with a velocity dispersion of cr = 278 km s ^ 
is shown in red. 



ent components of this complex system. In particular, while image 
iml may appear to be bluer than im2 from the magnitudes listed in 
Table [T] we suspect that this is just due to blending with the fore- 
ground red galaxy LI, because our X-shooter spectra of all three 
images, iml, im2, and im3, exhibit the same continuum slope. 

This paper is organised as follows. We describe the X-shooter 
observations in Section|2]and consider the lensing galaxies in Sec- 
tion[3] The main focus of the paper is the rich absorption and emis- 
sion line spectrum of the source, described in Sections |4] and |5] re- 
spectively. We discuss our main results in Section|6]and conclude in 
Section|7] Unless otherwise stated, we assume a flat cosmological 
model with Ho ~ 73 km s~^ Mpc~^ and = 0.72 jSpergel et al.l 
l2007b . 



2 OBSERVATIONS AND DATA REDUCTION 

The first of the secon d generation VUT instruments, X-shooter 
dD'Odorico et al. Il200d) is a three arm, single object echelle spec- 
trograph which was successfully integrated into standard VLT op- 
eration on October 1, 2009. The instrument covers simultane- 
ously the wavelength range from 300 nm to 2400 nm at resolv- 
ing powers R = A/AA = 5100, 8800 and 5600 in the UV-B 
(AA = 300-550 nm), VIS-R (AA = 550-1015 nm), and MR arms 
( AA — 1025-2400 nm) respectively. The corresponding slit widths 
are 1.0, 0.9, and 0.9 arcsec respectively. Although X-shooter is op- 
timised for observations of single targets with all arcsec long slit, 
it also allows for simultaneous recording of close pairs of objects 
on the sky. As a test of the performance of the spectrograph, two 
CASSOWARY candidates were observed during the commission- 
ing run of the instrument in March 2009. Result s on the firs t objec t, 
CSWA20, have recently been published by Pettin i et all ( |2010|) . 
This second paper reports the observations of CSWA5. 

Referring to Figure [T] the entrance slit of X-shooter was ro- 
tated to a position angle on the sky of P.A. = 8 degrees and aligned 
so as to capture the light of images iml and im2; at this PA, part 
of the light from image im3 also falls into the slit. The UV-B and 
VIS-R detectors were binned by a factor of two in the spectral di- 
rection. The total exposure time on source was 2400 s, split into 
two 1200 s long integrations; each integration on source was fol- 
lowed by a 1200 s exposure on a nearby blank sky region. From 
the r'-band acquisition frame reproduced in Figure[T] we measured 



a seeing full width at half maximum FWHM= 0'.'9. The seeing 
improved during the course of the observations; from the spatial 
extent of the spectra themselves we measured FWHM = 0'.'8 in the 
UV-B arm, 0'.'7 in the VIS-R arm, and 0'.'6 in the NIR arm. With 
the chosen setup, the slit losses in the three arms are estimated to 
be 15-10%. 

The spectra were reduced with a preliminary version of the 
ESO X-shooter pipeline ( Goldoni et al..2006,) , which uses subtrac- 
tion of the sky lines based on the procedures developed bv lKelsonl 
( |2003|) . For the UV-B and VIS-R data, we found it most advan- 
tageous to subtract the sky signal from the same two-dimensional 
(2D) frames on which the light from CSWA 5 was recorded, while 
for the NIR data sky subtraction from the separate blank sky 2D 
frames worked best. The pipeline reduction used calibration spectra 
taken during the commissioning run for order location and tracing, 
flat fielding and wavelength calibration. The final product from the 
pipeline is an extracted 2D, wavelength calibrated, rectified spec- 
trum with orders combined using a weighting scheme. For further 
data processing and analysis we used standard IRAF tools. We pro- 
duced ID spectra using a 1.2 arcsec extraction aperture in all three 
arms and for all three images of the source (see right panel of Fig- 
ure [TJ. As discussed in Section |4l the 16 strongest emission lines 
which were clearly detected in all three spectra were measured to 
be at the same redshift within Az — 5 x 10~^. Furthermore, we 
found no significant differences in the emission line ratios nor in 
the continuum slopes of the spectra of images iml, im2 and im3, 
leading us to conclude that all three images are lensed counterparts 
of the same background galaxy. Thus, unless otherwise specified, 
the following analysis was performed on the sum of the spectra of 
iml, im2 and im3, so as to maximise the S/N ratio of the data. 

Absolute flux calibration was performed with reference to the 
HST white dwarf standard GD 7 1 , while the 08 V star Hipparcos 
69892 provided a reference smooth spectrum for dividing out tel- 
luric absorption. Both stars were observed the night prior to the ob- 
servations of CSWA 5. Even though Hipparcos 69892 was observed 
at similar airmass as CSWA 5, the correction for atmospheric ab- 
sorption is not perfect, especially in the regions of low transmis- 
sion between the J, H, and K bands. A correction for atmospheric 
extinction was applied assuming a standard atmospheric extinc- 
tion curve. Overall, we expect the absolute flux calibration of these 
commissioning data to be accurate to within ~ 20%, based on the 
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comparison between the spectrum and the SDSS magnitudes of iml 
(the least contaminated of the four images). 



3 FOREGROUND GALAXIES 

Our X-shooter observations described above did not include either 
of the (presumably) foreground red galaxies visible in Figure [T] 
Galaxy L2 has an SDSS spectrum reproduced in Figure|2] We iden- 
tify a number of stellar absorption features, labelled in Figure |2l 
which indicate that it lies at redshift ^abs = 0.3877. In order to 
obtain an estimate of the velocity dispersion (and hence the mass) 
of the galaxy, we us ed the the penalised pixel fitting method of 
ICappellari & Emsellem ( 2004), taking into account the difference 
in velocity resolution be t ween the input stellar template spectra of 
ISanchez-Blazquez et alj ( |2006|) and the SDSS spectrum. The best 
fit gives a stellar velocity dispersion crL2 = 278 ± 6km s~^, 
where the uncertainty reflects the dispersion in the values of cr 
obtained with different stellar templates. The corresponding mass 
is M(r) = -Ko^r/G (appropriate to an isothermal sphere), or 
Ml2 2i 6 X 1O"M0, adopting r = lOkpc . 

At redshift z = 0.3877, the SDSS z filter is closest to rest- 
frame r. From the SDSS measured magnitude z{h2) = 17.67, 
we deduce an absolute magnitude A/, (L2) = —23.17, given the 
distance modulus m ~ AI = 41.17 in our cosmology and ap- 
plying an evolutionary correction of 0.852 to Afr, as proposed by 
[Sern ardi et al. (2006). There is no SDSS spectrum for galaxy LI. If 
it is at the same redshift as L2, z = 0.3877, its absolute magnitude 
is Mr (LI) = —22.82 and the projected separation between the two 
galaxies is 28kpc. The absolute magnitude and velocity dispersion 
of L2 are within the range of values found by Bernar di et al.l f200^ 
in their sample of massive early-type galaxies selected from the 
SDSS database. 

3.1 Lens Model 

For the interpretation of many of the properties of the source, it 
would be advantageous to obtain estimates of the magnification 
factors which apply to images iml-im3. To this end, we attempted 
a simple mo delling of the CS W AS system following the approach 
described in lPettini et al.l ( l2010h , but quickly realised that we have 
insufficient data at our disposal for even the simplest characteri- 
zation of the lens. For the idealised case of an isolated singular 
isothermal sphere, the typical deflection is given by: 

V200kms-V \Ds J 

where Dds is the angular diameter distance between deflector and 
source, whilst Ds is the distance between observer and source. If 
galaxy LI is: (a) the lens, and (b) at the same redshift as L2, then 
applying the measured aL2 = 278 km s~^ (which would actually 
be an upper limit for LI, since LI is fainter than L2 — see Table [B 
to eq. (|TJ would result in typical deflections of ~ 1.4arcsec, or 
only about half the deflections of images iml-im4. We are thus led 
to one of two conclusions. Either LI is substantially more massive 
than L2, as would be case if it is at a higher redshift, or L2 and 
other galaxies in the group significantly boost the lensing by LI. 
Presumably the second option is the more likely, given the proxim- 
ity of LI and L2 on the sky and their broadly similar magnitudes 

and colours (Table [TJ. 

Lensing by binary galaxies has been studied bv lShin & Evanj 

( l2008l) . From these authors' work, we find that three image mor- 
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Figure 3. Profiles of the [O II] AA3728, 3730 (upper panel) and 
[O III] A5008 (lower panel) emission lines extracted separately from im- 
ages iml, im2, and im 3. The spectra have been offset vertically, and scaled 
to facilitate comparison. The line profiles are indistinguishable within the 
accuracy of the data, indicating that the three images are lensed counter- 
parts of the same background galaxy at ^cm = 1.06856. 

phologies qualitatively corresponding to iml, im2 and im3 can be 
obtained (see lower panel of their Figure 10). However, five im- 
age morphologies qualitatively corresponding to iml, im2, im3 and 
im4, together with a highly demagnified fifth or central image are 
also possible (see middle panel of their Figure II). Thus, in order to 
make progress, it will be necessary to measure the redshifts of LI 
and other galaxies in the field, as well as to resolve the ambiguous 
status of im4. Without a satisfactory lensing model, any luminosity- 
dependent quantity derived for the source will be uncertain by an 
amount corresponding to the unknown magnification factor. 



4 SOURCE GALAXY 

Our X-shooter observations show images iml, im2 and im3 of the 
CSWA5 system to be gravitationally lensed images of a luminous 
star-forming galaxy at z = 1.0686. This is illustrated in Figure[3] 
where the individual profiles of two of the best observed emission 
lines, the [O II] AA3728, 3730 doublet and the [O III] A5008 line 
are shown separately for images iml, im2 and im3. The three sets of 
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Figure 4. Examples of interstellar absorption lines in the X-shooter spectrum of CSWA 5. The full list of lines identified is given in Table[2] 



Table 2. Interstellar absorption lines in the X-shooter spectrum of CSWA 5. 



^obs (A) 


Ion 


^'ab (A) 


2;* 
abs 


Wo (A) 


log(Af/cm-2) 


3740.44 


Sill 


1808.013 


1.06881 


0.47±0.20 


14.7 


3835.81 


Aim 


1854.520 


1.06836 


0.83±0.27 


>12.9 


3853.27 


Aim 


1862.790 


1.06855 


0.61±0.20 


>13.2 


4848.81 


Fell 


2344.214 


1.06842 


1.54±0.18 


>13.8 


4912.11 


Fell 


2374.461 


1.06873 


1.07±0.11 


14.0 


4928.75 


Fell 


2382.765 


1.06850 


1.14±0.09 


>13.0 


5349.78 


Fell 


2586.650 


1.06823 


0.94±0.15 


>13.7 


5378.10 


Fell 


2600.173 


1.06836 


1.32±0.16 


>13.3 


5783.56 


Mgll 


2796.351 


1.06825 


2.87±0.6 


>13.4 


5798.96 


Mgll 


2803.531 


1.06845 


3.06±0.7 


>13.7 


5901.91 


Mgl 


2852.964 


1.06869 


1.23±0.6 


>12.4 



Vacuum heliocentric. 
^ Vacuum rest wavelengths. 



profiles are indistinguishable within the noise. Furthermore, when 
measured separately in each image, the redshifts of the 16 strongest 
emission lines agree to within Az = 5 x 10~^. On the basis of this 
evidence, we conclude that iml, im2 and im3 are gravitationally 
lensed images of the same galaxy. 

The galaxy exhibits a rich absorption and emission line spec- 
trum which we now describe in detail. Unless otherwise specified, 
the measurements reported below refer to the sum of the three im- 
ages, so as to maximise the S/N ratio. 



4.1 Absorption Lines 

The UV-B and VIS-R spectra of the source show many strong ab- 
sorption lines from neutral and singly ionised species formed in the 
interstellar medium of the lensed galaxy. They are listed in Table[2l 
and a few examples are reproduced in Figure|4l They define a mean 
absorption redshift Zabs ~ 1.0685±0. 0002, where the error quoted 
is the standard deviation from the mean. The absorption lines de- 
tected are wide, spanning a velocity range of up to ~ 400 km s^^. 

With the relatively short exposure time devoted to CSWA 5 
during the X-shooter commissioning, the signal-to-noise ratio of 
the spectra is low in the continuum, S/N~ 5. Thus, only the 



strongest interstellar absorption lines could be detected with con- 
fidence. Such lines are saturated and unsuitable for column density 
determination. Under these circumstances, application of the appar- 
ent optical depth method o f ISavage & Sembachl i ll99l]) yields lower 
limits to the column densities A'^. Out of the transitions listed in Ta- 
ble [2l only two, SillA1898 and FeIlA2374, may be sufficiently 
weak to provide estimates of A'^(Sill) and A'^(Fell) respectively. 
Taken at face value, the column densities of these two ions may 
indicate a depletion of Fe by a factor of ~ 5 , since Si and Fe hav e 
comparable relative abundances in the Sun jAsDlund et Zll2009h . 
whereas N{¥e II) ~ 1 /57V(Si II) in CSWA 5. 



4.2 Emission Lines 

Thanks to its wide wavelength coverage, the X-shooter spectrum 
of CSWA 5 includes a multitude of nebular emission lines, from 
He I A10833 to the C III] AA1907, 1909 doublet; a selection is re- 
produced in Figure |5] As can be seen from Table [3] we detect 
11 members of the Balmer series, as well as Pa(5, and several 
transitions of He, C, N O, Ne, S and Ar in a variety of ioni- 
sation stages. The emission lines define a mean redshift Zcm = 
1.06856 ± 0.00004 (la), weighted by the flux in each line. 

The strongest emission lines are best fit with a combina- 
tion of two Gaussian functions; a narrow component centred at 
zi = 1.06859 with a velocity dispersion ai = 50 ± 4km s~^ 
(after correction for the instrumental resolution) and a second, 
broader component centred at the lower redshift Z2 = 1.06843 
with (72 = 120 ± 10 km s^^. The uncertainties reflect the spread 
of values from acceptable fits to the line profiles. The redshift dif- 
ference corresponds to a velocity separation Av — 23 km s~^. The 
top left-hand panel in Figure [5] illustrates the profile decomposi- 
tion for the [O II] AA3728, 3730 doublet. The flux ratio between the 
narrow and broad component varies between 1.5:1 and 2:1 among 
different emission lines. The two emission components are not eas- 
ily distinguishable in the weaker, noisier lines. Thus, in Table[3]we 
quote the integrated fluxes across the emission lines, which are in- 
dependent of profile decomposition. 
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Figure 5. Examples of nebular emission lines in the X-shooter spectrum of CSWA5. The top left-hand panel illustrates the profile decomposition into 
two components for the [O II] AA3728, 3730 doublet. Black histogram: observed profile of the [O II] emission lines. Red dotted lines: individual Gaussian 
components; each member of the doublet is fitted with two emission components with the parameters given in Section ITzl Red continuous line: total fitted 
profile to the doublet. The blue histogram near zero intensity shows the difference between the observed and fitted profiles. 



5 PHYSICAL PROPERTIES 

The emission line fluxes collected in Table [3] and their ratios al- 
low us to deduce a number of physical parameters characterising 
the H II regions of the lensed galaxy in CSWA5. In the following, 
we make the assumption that the gas is ionised by OB stars, with 
negligible contributions from shocks and an active galactic nucleus 
(AGN; see SectionlOl. 

Physical quantities that are deduced from the ratios of emis- 
sion lines — such as reddening, gas temperature, density, and chemi- 
cal composition, or from the relative spectral energy distribution — 
such as the age of the stellar population, are independent of the 
magnification factor applicable to each of the three images covered 
by our observations. This is not the case, however, for quantities 
derived from line and continuum luminosities, such as the star for- 
mation rate (SFR) and stellar mass (M*). In placing these quanti- 
ties on an absolute scale, we are hampered not only by the lack of 
a lensing model (Section[3TTJ, but also but the unknown fraction of 
the light from image im3 captured by the X-shooter slit (see Fig- 
ure [TJ. To simplify matters, we used the spectrum of image iml 
alone, rather than the sum of iml, im2 and im3, in our estimates of 
the star formation rate (Section r5.2l ( and stellar mass (Section[53](. 
We chose iml because it is less affected by blending than the other 
images (see FigurefTJ, and we include explicitly the unknown mag- 
nification factor l//icns in our estimates of SFR and A'h below. 
However, we suspect that /icns may be of order unity (within a fac- 
tor of 2-3) for a single image, as found by Pettini et al. (2010) for 
CSWA20, which has a similar configuration to CSWA5. 



5.1 Reddening 



The observed ratios of Balmer line intensities can be compared 
with theoretical expectations from recombination theory to deduce 
the amount of dust reddening. As is common practice in the anal- 
ysis of H II regions, we assumed Case B recombination, an elec- 
tron temperature Te = 1 0^ K, and el ectron densities in the range 
n(e) = 10^-10'' cm"^ tosterbrockl 1 1989.) . We did not include 
Hq in the analysis because at z = 1.06856 it is redshifted to 
Aobs ~ 13 579 A, between the infrared J and H bands, where 
atmospheric absorption is very severe and difficult to correct for 
accurately (which presumably explains why the observed ratio 
Ha/H/3 in Table[3]is ~ 25% lower than the value Ha/H/3 = 2.86 
expected for Case B recombination). From the observed H/3/H7 ra- 
tio, we deduce E{B - V) = O.W ± 0.03, while the H/3/H(5 ratio 
gives E{B — V) = —0.09 ± 0.05. The higher order Balmer lines 
have larger relative line flux uncertainties, so the reddening derived 
from higher order ratios is less certain. The weighted average is 
E{B -V) = 0.03 ± 0.02. Thus, the Hll regions of CSWA5 (or 
at least the unobscured portions we see) suffer little dust redden- 
ing. Since the extinction is low (A\r ~ 0.1 mag, taking the usual 
ratio of total to sel ective extinction Ry = Av/E{B — V) ~ 4, 
ICalzetti et aLll2000l) and consistent with zero within the uncertain- 
ties, we did not include it in the analysis described in the following 
subsections. 
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Table 3. Nebular emission lines in the X-shooter spectrum of CSWA5. 
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434]^ 69 
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10059.59 


Hl(H/3) 


4862.69 
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10261.48 


[OlII] 


4960.30 


1.06872 


16.84±0.27 


10369.68 


[Olll] 


5008.24 


1.06873 


49.53±0.45 


12159.93 


He I 


5877.25 


1.06898 


1.20±0.15 


13037.30 


Ol 
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1.06874 
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H I (Ho) 


6564.66 
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1.06877 


4.33±1.13 


14618.95 


He I 


7067.14 


1.06858 


0.12±0.05 


14765.79 


[Arm] 


7137.77 


1.06868 


0.17±0.05 


19722.90 


[SIII] 


9533.71 
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10052.2 


1.06947 
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He I 


10833.3 
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^ Vacuum heliocentric. The redshift errors are typically 2 X 10~^, 
although it the worst cases — where a feature is contaminated by sky 
residuals or telluric absorption, the redshift error can be up to ten times 
higher. 

^ Vacuum rest wavelengths. 

Integrated line fluxes in units of 10~^^ erg s~^ cm~^. The errors 
quoted reflect only the random uncertainties from the counting statistics 
and do not include systematic sources of eiTor resulting from connection 
for atmospheric absorption, subtraction of sky emission lines, and 
absolute flux caHbration. For some emission lines (e.g. Ha), the 
systematic errors can be much larger than the random errors. 



5.2 Star Formation Rate 

We can obtain estimates of the star formation rate from the lumi- 
nosities in the Ha emission line and t he near-UV contin uum, using 
the calibrations of these measures bv lKenmcut3 ( ll998h . Using the 
H/3 flux measured for iml, F(H^) = 3.1 x 10"^^ erg s~^ cm~^ 
and the unreddened Case B recombination ratio F(HQf) /F(H/3) = 



2.86, we deduce F(Ha) = 8.9 x 10"^'' erg s"^ cm"^ and a lu- 



minosity L(Ha) 



5.1 X 10*^ erg s ^. This in turn implies: 



SFR = 7.9 X 10~*^ L{Ha 



1 1 
X — X - — 

1-8 /ion; 



23(M0 yr-^), (2) 



using the conv ersion between L(Ha) and SFR proposed by 
iKennicutj ( ll99S) . The two corrections factors account for: (i) the 



flatten ing of the initial mass function (IMF) b elow 1 M p) I Chabriej 
l2003h relative to the single power law of the Salpete 3 il955) IMF 
assumed bylKennicutt (1998), and (ii) the unknown lensing magni- 
fication of image iml, /ions, as discussed above. 

Turning to the UV continuum, we measure F — 6.1 x 
10~^^ erg cm~^ s~^ from our X-shooter spectrum of iml at 
observed wavelengths near 5800 A. The corresponding rest-frame 



luminosity near 2800 A, Li,(2800) = 1.9 x 10^^ erg s"^ Hz"\ 
implies: 

SFR= 1.4x10"^* L^(2800)x-i-x--^ ~ 15(M0 yr"^),(3) 

1-8 /ions 

with the same correction factors as above. The two values of SFR 
are in reasonable agreement, given the systematic uncertainties af- 
fecting the two estimators, which sample different portions of the 
IMF (see .Meurer et al..2009,) . 



5.3 Electron Temperature and Density 

At the S/N ratio of the present data, the auroral lines most 
commonly used for the determination of the electron temper- 
ature, [0III]A4363, [Nll]A5755, and [SlIl]A6312 are below 
the detection limit. The 3a upper limit F(4363) ^ 0.3 x 
10"^'' erg s~^ cm~^ tra nslate s to an upper limit on the temper- 
ature To < 9900 K jAllej Il984l) . Turning to density sensi- 
tive line ratios, we resolve both the [O II] AA3728, 3730 and the 
CIII]AA1907, 1909 doublets, while [S II] AA6716, 6731 falls be- 
tween the J and the H bands, where the atmospheric transmission 
is a few percent only. We used the NEBULAR package in IRAF, in- 
corporating the calculations bv lShaw & Dufourl ( 1 19951) . to deduce 
n(e) = 315 ± 15 cm"^ from the observed F(3726)/-F(3728) ra- 
tio, assuming To = 10 000 ± 1000 K. Although less tightly con- 
strained, the observed F(1907)/F(1909) ratio is consistent with 
the electron density deduced from the higher S/N ratio [O II] dou- 
blet. 



5.4 Chemical Abundances 

With many nebular lines detected, we can estimate the abundances 
of the C, N, O group using a variety of strong emission line diag- 
nostics. Results are collected in Table|4] 

Oxygen. Oxygen is the element most commonly used to char- 
acterise the overall degree of metal enrichment of Hll regions 
and we thus consider it first. Among the many calibrations of the 
oxyg en abundance based on the ratios of strong emission lines 
(e.g. iKewlev & Ellisonll2008l) . the 7?23 method first proposed by 
|Pageletal.Ul979l) is the one still most commonly applied, at least 
to low redshift galaxies. We used the calibration of the index 
R 23 = [F(3728) + F(3730) + F(4960) + F(5008)]/F(H/3) 
bv lKobulnic kv & Ke wlevI ( |2004|) . appropriate to the lower branch 
of the double-valued J?23 vs. 0/H relation, to deduce 12 + 
log (0/H)^23 — 8.28, or approximatel y 40% of the oxygen abun- 
dan ce in the Sun I Asplund et al. I I2OO9 I) and in the Orion nebula 
( Este ban et al The random error in the value of log (O /H) 

is small and negligible compared with the 0.2- 0.3 dex systematic 
uncertainty in the calibration of the 7?23 index iKewlev & EUiso j 
I2OO8I) . The lower branch solution is favoured by the low value of 
the N abundance deduced below. 

The calibrations of the oxygen abundance bas ed on the ratios 
of nitr ogen, oxygen and Balmer lines proposed bv lPettini & Pagell 
1 I2OO4I) have the advantages, compared to R23, of being single- 
valued and relatively robust to errors in flux calibration and 
dust extinction. However, in the case of CSWA5, their use is 
possibly compromised by the awkward redshifted wavelengths 
of Ha and the [N II] AA6550, 6585 doublet which fafl between 
the J and H bands, as already explained. We deduce 12 + 
- 8.47 ± 0.07 from consideration of the ratio 



log (O/H), 



N2 



\og[F(6585)/F{Ha)], and 12 + log(0/H)^ 
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Figure 6. X-shooter spectrum of image iml in CSWA5, colour-coded to show tlie spectral ranges covered by the three arms of the spectrograph, as follows: 
blue is UV-B, green is VIS-R, and red is NIR. The spectrum has been smoothed with a boxcar filter of width 2.5 A. Overplotted in black is the best- 
fitting template spectrum computed with a lBruzual & Charloll i2003l) model for a 40 Myr-old continuous star formation episode and assembled stellar mass 
M* = 2.8 X IO^Mq. 



Table 4. C, N, O Abundances in CSWA5. 



12 + log{0/H) 


log(N/0) 


log(C/0) 


8.3-8.5 


-1.0 ±0.1 


-0.6 ±0.2 



8.29 ± 0.04 from the ratio 03iV2 = log[F(5008)/F(_ff/3)] - 
log[F(6585) /F(H a)\, using in both cases the calibrations of these 
indices deduced byl Pettini& Pagell ( |2004|) . A gain, the random er- 
rors quoted are small compared to the 0.2-0.3 dex systematic un- 
certainties of the calibrations. We conclude that the oxygen abun- 
dance in the lensed galaxy in CSWA5 is 12 + log (0/H) = 8.3- 
8.5, or (0/H)cswA5 ^ 0.4-0.65 (O/H)^. 

Nitrogen. In the absence of temperature- sensitive aurora l lines , 
we make use of the method developed by iThurston et al.l j 19961) 
to deduce the relative abundances of N and O. These authors used 
photoionization models to deduce a relationship between the [N II] 
temperature and the i?23 index; in our case, f[]vii] = 11000 ± 



1000 K. Once the [N ll] temperature is known, the N/O ratio can 
be deduced directly from the r atio of the [N ll] AA6550, 6585 
and [O II] AA3728, 3730 doublets jPagel et al.lll992h. under the as- 
sumption that N+/0+ ~N/0 which IThurston et al.l ( Il996h argue 
introduces only a small error. Following this procedure, we find 
log (N/0)cg^A5 = - 0.95 ± 0.10 or ~ 80 % of the solar value 
log(N/0)g = -0.86 jAsplund et a"i]|2009h . A slightly subsolar 
N/O ratio when the abundance of oxygen is ~ 0.5 solar is in line 
with the large body of such measuremeiits novy available in galax ies 
at a range of redshifts (e.g. llzotov et al.ll2006l;|Pettini et alj2008l) . 

Carbon. In the low density limit, the ratio C'^'^IO'^'^ can be 
deduced from the measured fluxes of the C III] AA1907, 1909 
and [O III] AA4960, 5008 doublets and a knowledge of the rel- 
evant collision strengths. Using the line emissivities given by 
Aller (1984), we deduce log(C+VO+^) = -0.6 ± 0.2, for 
Tc = 10 000 ± 1000 K. The ionization correc tion factor to deduc e 
C/0 from C+^/0+^ is expected to be small jOarnett et al]|l995[) . 
and therefore we find ( C/0)csvyA5 — 0.3 -0.7(C/O)(;), where 
log(C/O)0 = -0.26 jAsplund et a"i]|2009l) . As for nitrogen, the 
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sub-solar C/0 ratio in CSWA5 conforms to the established be- 
havio i ir of C /0 vs. O/H in metal-poor H II regions (Garnett et al. 



I995L ll999S. Galactic stars /Akerman et al.' 2004; Fabbian et al. 



2010) 



200S ), and Lyman break galaxies dShaplev et al.. .2003. : .Erb et al.. 



5.5 Age and Assembled Stellar Mass 

Finally, we can obtain estimates of the age and stellar mass of the 
lensed galaxy by comparing its observed spectral energy distribu- 
tion (SED), from the rest frame UV to the near-IR (from ~ 2000 A 
to ~ 1.2/^m where our X-shooter spectrum has the highest S/N ra- 
tio), to those of synthetic spectra computed with simple stellar pop- 
ulation models (e.g. Erb et al. 20 06b). For this purpose, we used 
the popu l ation synthesis code of Isruzual & Charlo t| ( I2OO3I) with 
IChabrierl j2003l) IMF, metallicity Z = 0.4^0 (Section [531l, and 
no internal reddening (Section f5.lt . We generated two families of 
models with these parameters, for the two limiting cases of continu- 
ous star formation and an instantaneous burst, in each case varying 
the age from 1 to 70Myr. The best fitting models and the corre- 
sponding values of age and assembled stellar mass were determined 
by minimizing the value of x^, given by; 



(/o 



X /model, A ) 



(4) 



where fobs.x and (Tobs.A are the observed flux at wavelength 
A and its erro r respectively, and /model, a is the flux of the 
iBruzual & Ch ariot ( 2003) model spectrum at the same wavelength. 
The normalization factor b gives the stellar mass of the galaxy. 

In fitting the model spectra to the X-shooter spectrum of image 
iml, we excluded regions affected by strong telluric absorption or 
by prominent residuals in the subtraction of sky emission lines. We 
found that the best fitting models (see Figure[6} are those with ages 
of 40 ± lOMyr and stellar masses M. = (2.8 ± 1.0) x 10^ x 
l//icns M0, where /icns is the unknown magnification factor of 
iml, and the error includes the 20% uncertainty in the absolute flux 
calibration (Section[2l(. 

These values appear to be robust to the choice of star forma- 
tion mode, with the single burst and continuous star formation mod- 
els converging to similar solutions. We also investigated the possi- 
bility that the stellar continuum suffers a greater extinction than 
the Hll emission lines, since the former could in principle sam- 
ple a different (presumably older) stellar population than the latter. 
However, without recourse to the rest frame far-UV spectral range 
(which is inaccessible from the ground at z = 1.0686), we come 
up against the well-known age-extinction degeneracy. For example, 
models in which the stellar continuum is reddened with a colour ex- 
cess E{B — V) — 0.3 also provide satisfactory fits to the X-shooter 
spectrum of CSWA5, albeit with younger preferred ages of 25 ± 
7 Myr and lower stellar masses of (2.3 ±0.7) X 10^ X l//lens Mq. 

We can obtain an independent estimate of the age of the star- 
burst from consideration of the equivalent width of the H/3 emission 
line, WoiUP) ~ 100 A, using the values of WoiUp) as a func- 
tion o f time calculated wit h the STARBURST99 spectral synthesis 
code jLeitherer et al.ll2001 ). The equivalent widths of the B aimer 
recombination lines fall rapidly following a burst of star formation 
and Wo (H/?) < 100 A for all ages greater than 5 Myr. On the other 
hand, the time dependence is less steep for continuous star forma- 
tion and ages of ^ 25 Myr are indicated, in better agreement the 
age inferred from the analysis of the UV to near-IR SED. 



6 DISCUSSION 

Summarizing our findings, we have established that CSWA 5 is an 
actively star- forming galaxy at z = 1.0686 lensed by an apparent 
foreground group of massive red galaxies, at least one of which 
is at z = 0.3877. The lensed source is forming stars at a rate 
SFR ~ 20 Mq yr~^, suffers negligible reddening, has an oxy- 
gen abundance of approximately half-solar, and sub-solar N/0 and 
C/O ratios (as expected for low metallicity galaxies). Its blue spec- 
tral energy distribution, from the rest-frame UV to the near-IR, and 
high H/3 equivalent width are indicative of a young age, only ~ 25- 
50 Myr, during which time the galaxy has assembled a stellar mass 
of '--^ 3 X W^Mq. The star-formation rate and stellar mass are un- 
certain by an unknown magnification factor, which may be of order 
unity for each of the four gravitationally lensed images of CWSA 5. 

These physical characteristics are broadly in line with those of 
the population of UV-selected galaxies at redshifts z — 1-3 (e.g. 
IPettini et al.ll2007l . and references therein), as we now discuss. 



6.1 Mass-Metallicity Relation 

The most relevant compari son here i s with the mass-metallicity re- 
lation determined bv Savaglio et al. I ( I2OOS) for galaxies at a mean 
(2:) ~ 0.7 selected from the Gemini Deep Deep Survey (GDDS) 
and shown in Figure [7] It must be borne in mind here that such 
comparisons are fraught with pitfalls for the unwary, because of 
systematic differences between the methods used to determine both 
the stellar mass and the m etallicity. In this particular case, the study 
by Savag lio et al. should be compatible with ours, as both 

used a lChabrieJ i2003h IMF in arriving at the stellar masses and the 
R23 index for the oxygen abundance (although the reliability of 
the latter at app arently super-solar metallicities is questionable — 
see |Pettinil2"00a) . 

An additional complication is the unknown magnification fac- 
tor which applies to our derivation of the stellar mass. If /lens in 
Section 15.51 is of order unity, then Figure [7] shows that CSWA 5 
is somewhat metal-poor for its stellar mass, compared to GDDS 
galaxies at z ~ 0.7. Definite conclusions are difficult at this stage, 
given the scatter of the GDDS galaxies about the mean relation in 
Figure [T] and the lack of objects as m etal-poor as CSWA 5 in the 
sample considered by I Savaglio et alj l|2005,). Taken at face value, 
the offset of CSWA 5 in Figure |7] is in the same sense as the more 
ge neral redshift evolution of the mass-metallicity relation proposed 
bv lMaiolinoetal.l ( l2008h . 



6.2 Specific Star Formation Rate 

Since the lens magnification factor enters into the calculation of 
both SFR and M,, the ratio of these two quantities, commonly re- 
ferred to as the specific star formation rate (SSFR), should not be 
affected. For CSWA 5 we deduce a specific SFR = (8.2 ± 2.2) x 
10~^ yr~^. The reciprocal of this value gives a timescale for the 
star formation activity of ~ 120 ± 45 Myr, which is ~ 3 times 
higher than the age deduced in Section [53] from SED fitting. 

In Figure [8] we compare the value of SSFR in CS WA 5 with 
those of UV-selected galaxies at 2 ~ 2 from the work bv lErb et"ail 
(2006^ and of lower redshift galaxies from the VVDS survey by 
iLamareille et alj ^2W% . Clearly, CSWA 5 is more akin to the for- 
mer than to the latter in its relatively high SSFR and correspond- 
ingly short star formation timescale. 
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Figure 7 . Mass metallicity rel ation for galaxies at 2 ^ 0.7 from the GDDS 
fcrosses JSavaglio et al.ll2003) . The location of CSWA5 is indicated by the 
large filled circle. The errors shown do not include the systematic uncer- 
tainties in the metallicity and mass determinations. 



6.3 Emission line ratios 

Figure [5] shows that CSWA5 is well separated from the lo- 
cus of SPSS galaxi es in the ionization diagnostic diagram of 



[Baldwin et al.l 1 198 ih . sometimes referred to as the BPT diagram. 



While we find no evidence for the presence of an AGN in CSWA5 
from the width of the emission lines, we note that such offsets now 
app ear to be common for the most a c tively star-forming galax- 
ies terb et alJ l2006al : iLiu et alJ l2008h . iBrinchmann et all d2008l) 
pointed out that, among SDSS galaxies, there is a correlation be- 
tween excess SSFR and offset from the main locus occupied by 
'normal' galaxies in the BPT diagram. Galaxies with higher SSFRs 
tend to exhibit higher values of the ratios [O ni]/H/? and [N II]/Hq. 
CSWA5 apparently fits this trend. The underlying physical expla- 
nation for the offset is not well established, but there are indications 
that it may be related to a higher ionization parameter of the H II re- 
gions, w hich in turn is p r esum ably linked to the high star formation 
activity jLiu et all2008l ; lBrinchmann et al.ll2008h . 



7 CONCLUSIONS 

It is encouraging and satisfying that through dedicated searches of 
large area surveys the number of gravitationally lensed galaxies at 
high redshifts continues to increase. The work presented here has 
not only identified a new such example in the CASSOWARY cat- 
alogue but, more importantly, has demonstrated the wealth of in- 
formation on the physical conditions of these objects which can be 
gleaned from even relatively short exposures with the X-shooter 
spectrograph, thanks to its high efficiency and wide wavelength 
coverage. The main impediment we have found to the full char- 
acterization of CSWA 5 is the lack of a tensing model which can 
reproduce the deflections of the images. Higher spatial resolution 
imaging of the field (achievable with the Hubble Space Telescope) 
and more redshift determinations of galaxies in the field are needed 
for a comprehensive census of the mass distribution along the line 
of sight and an accurate estimate of the magnification of CSWA 5. 
Such information, combined with higher S/N ratio spectra of the 
four lensed images (the work presented in this paper is based on 
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Figure 8. Specific star-formation rate of CSWA 5 (large black dot) com- 
pared with those of galaxies at z > 2 (grey dots; lErb et al]|2006 b). and 
galaxies at < z < 0.6 from th e VIMOS VLT Deep Survey (VVDS) 
(small dots; iLamareille et aljr2009h . All SFRs were derived from the Ha 
luminosities. Lines of constant SFR are indicated by the straight solid and 
dashed lines. The mass (but not the specific SFR) of CSWA 5 is uncertain 
by an unknown factor 1 / /icns which may be of order unity. The scale on 
the right vertical axis shows the formation time scale (the inverse of the 
specific SFR). 



only 40 minutes of integration on target), would lead to a full map- 
ping of the properties of this star-forming galaxy. 

As a final note, we should like to stress that the use of gravita- 
tional lensing to access the properties of galaxies at high redshifts 
is in some way still in its infancy. While the few cases studied so 
far have been very successful in providing detailed views of these 
systems which would otherwise have simply been unobtainable, it 
is also the case that nearly all of them appear to be drawn from the 
bright end of the galaxy luminosity function. Fainter galaxies are 
accessible with 8-10 m telescopes equipped with modem spectro- 
graphs, but are often overlooked in searches for strong gravitational 
lenses. By systematically extending such searches to fainter limits 
with current and future (e.g. with the Large Synoptic Survey Tele- 
scope) surveys, it will be possible to identify intrinsically fainter 
galaxies, particularly at epochs when the Universe was forming 
most of its stars, and thereby assemble a full picture of the galaxy 
formation process in action. 
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